float f(float x){
  return x*x-5*x+1;
}